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ZERO-GRAVITY VENTING OF THREE REFRIGERANTS 
by Thomas L Labus, John C. Aydelott, and Geraldine E. Amling 

Lewis Research Center 


SUMMARY 

An experimental investigation of venting cylindrical containers partially filled with 
initially saturated liquids under zero-gravity conditions was conducted in the NASA 
Lewis Research Center 5-second zero-gravity facility. The test fluids possessed a 
near 0° contact angle on the cylindrical container surface resulting in a hemispherical 
liquid -vapor interface shape in zero gravity. The effect of interfacial mass transfer on 
the ullage pressure response during venting was analytically determined. This pressure 
response was compared with both the experimental pressure response and the pressure 
response based on an adiabatic decompression computation. The results showed that 
interfacial mass transfer, based on a conduction analysis applied to an infinitely planer 
(flat) liquid -vapor interface, in a lumped system analysis is significant in determining 
the ullage pressure response. The pressure response computations from adiabatic de- 
compression consistently yielded too large a predicted ullage pressure decrease. In- 
clusion of the interfacial mass transfer effects resulted in approximately a 30 -percent 
improvement in the predicted pressure response. 


INTRODUCTION 

The use of high-energy liquid propellants in our space program has led to a need 
for information concerning the thermodynamic behavior of cryogenic fluids in tanks 
which are vented or depressurized in space. The task of venting in low gravity has been 
successfully accomplished during a number of past missions with venting systems that 
rely exclusively on auxiliary thrusters to actively position the liquid propellant away 
from the tank vent. 

Short-term venting has been employed by the Saturn S-IV-B (ref. 1) and the Centaur 
(ref. 2) vehicles in near-Earth orbit. The methods of pressure control include using the 
vented propellant vapor to provide adequate acceleration to keep the liquid propellant 
settled at the end of the tank opposite the vent or using auxiliary chemical rockets to 
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provide a reorientation maneuver to ’’collect” the liquid propellant at the end of the tank 
away from the vent. The low-level acceleration criteria for controlling the liquids are 
based on studies conducted in the zero-gravity drop tower (ref. 3). 

The problems of long-term space missions such as interplanetary flights, deep space 
probes, and orbital storage systems more acutely point up the need for basic under- 
standing in the area of zero- and low -gravity venting. The short-term cryogenic storage 
requirements of the Saturn and Centaur vehicles yield optimum weight system designs 
when the cryogenic boiloff is sacrificed in lieu of additional insulation. For long-term 
storage of cryogenics, insulation systems will be improved such that the vented propel- 
lant vapor will not be sufficient to supply propellant settling, and the use of auxiliary 
thrusters will place too large a weight penalty on the vehicle. Unless a device such as 
a liquid propellant thermal conditioning system (ref. 4) is employed, propellant control 
becomes mandatory. However, the application of a thermal conditioning system may add 
excessive weight penalties and complexity. 

If some type of internal surface tension devices are included in the propellant tank 
to maintain liquid-vapor interface control, it would appear that venting during long-term 
space missions could be considered on a continuous basis. If the vent rate is very low, 
pressure control may be accomplished without disturbing the liquid bulk since liquid- 
vapor interface evaporation alone will supply the required vapor for venting. The ob- 
jective of this study is to predict the pressure response of a saturated liquid-vapor sys- 
tem when undergoing a venting or depressurization process in zero gravity (weightless- 
ness) at low vent rates. 

A venting analysis was formulated based on lumping the continuity and energy equa- 
tions governing the vapor space. The analysis is not limited to any particular tank shape 
or contained liquid. 

The analysis includes interfacial mass transfer based on an infinitely planer (flat 
surface) conduction analysis (refs. 5 and 6). Interfacial mass transfer is often neglected 
for rapid venting, but for low vent rates the effect of interface mass transfer on the 
pressure response is significant. The pressure response for rapid venting is usually 
based on the adiabatic decompression venting model. 

The pressure responses determined from the venting model which includes inter- 
face mass transfer are compared with the results from the adiabatic decompression 
venting model and the pressure responses obtained from short duration (5 sec nominally) 
drop tower tests conducted at the Lewis zero-gravity facility. Some qualitative informa- 
tion concerning bulk boiling similar to that obtained in reference 7 is also presented. 

The work presented herein is concerned primarily with low vent rates in contrast to 
reference 8 in which larger vent rates (more than one ullage/sec) were of interest. Ref- 
erence 8 presents experimental zero-gravity data for both refrigerant Freon 11 and 
liquid hydrogen. The test fluids used for this study include the refrigerants Freon 11 


2 



(Rll), octofluorocyclobutane (RC318), and n-butane (R600). These fluids were chosen 
because they provide a range of heat -transfer properties and are saturated at room tem- 
perature and pressures only slightly above atmospheric. 
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SYMBOLS 

o 

area, m 

discharge coefficient 

specific heat at constant volume, J/(kg)(K) 
functions 

specific enthalpy , J /kg 
heat of vaporization, J/kg 
thermal conductivity, W/(m)(K) 
mass, kg 

unit normal vector 

O 

pressure, N/m 

volumetric flow rate, m /sec 

heat flux , W / m^ 

ga s con stant , (m) (N) /(kg) (K) 

temperature, K 

time, sec 

internal energy, J 

specific internal energy, J/kg 

volume, m^ 

velocity, m/sec 

thermal diffusivity, m /sec 

ratio of specific heats 

density, kg/m 

parameter 
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Subscripts: 

1 initial state 

2 state for times, t > 0 

adia adiabatic decompression computation 

anal analytical including interface mass transfer 
exp experimental 

f final 

l liquid 

o vented vapor 

s liquid -vapor surface 

sat sa to at ion conditions 

v contained vapor 


APPARATUS AND PROCEDURE 


The Lewis zero-gravity facility was used to obtain the experimental data for this 
investigation. A complete description of the facility, the experiment package, and pro- 
cedures for conducting the tests can be found in appendix A. The venting tests were 
conducted with right -circular cylindrical containers having flat ends as shown in fig- 


ure 1. These test containers were fabricated from acrylic plastic. 

The liquids employed in the tests were refrigerant 11 (CCl^F), refrigerant C318 
' C 4 F 8' » and refrigerant 600 (CHgCHgCHgCBg) . Refrigerant 600 is commonly known as 
n-butane and refrigerant 11 is commonly known as Freon 11 or Genetron 11. The per- 
tinent thermodynamic properties for all these fluids, such as enthalpy, specific volume, 
and entropy, can be found in standard refrigeration tables (ref. 9). All three of these 
fluids were found to exhibit a nearly 0° static contact angle on the test container sur- 
face, and are, therefore representative of typical liquid propellants. The area of the 
liquid surface (5. 57X10 m ) is simply the area of a hemisphere having a radius equal 
to that of the tank since the venting sequence occurs during weightlessness. The tank 
pressure and the temperature in the vicinity of the liquid -vapor interface were moni- 
tored continuously before the drop test to ensure reaching nearly saturated conditions. 
The tank pressure was recorded during the drop test along with high-speed motion pic- 
tures. After releasing an experiment package, a maximum of 1. 9 seconds was allowed 
for the liquid to achieve the hemispherical zero-gravity interface configuration; the 
timers then activated the vent sequence for approximately 3 seconds. 
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ANALYSIS 


Adiabatic Decompression Venting Model 

The pressure reduction of a tank containing only a perfect gas during the adiabatic 
withdrawal of that gas is calculated in order to make comparisons with a venting model 
which includes interfacial mass transfer. Conservation of mass is applied in the follow- 
ing manner: The rate of decrease of mass in the tank is set equal to the withdrawal rate 
where the volumetric withdrawal rate Q is assumed constant. When proceeding along 
the lines of the analysis presented in appendix B, substituting equation (B16), which is 

P 1 T 2 
— - - e 


into equation (Bll), which is 


Tj _ e -QRt/V v C y 


leads to the following expression for the pressure reduction during an adiabatic decom- 
pression: 


^2_ e ‘ yQt/V v 


(1) 


Interface Mass Transfer Venting Model 

At time zero, the venting process is assumed to start with initially saturated vapor 
exiting through the vent to a vacuum. A schematic of the venting system at various in- 
stances of time is shown in figure 2. The initial thermodynamic state of the vapor and 
liquid are known (saturated conditions). The pressure decrease within the container dur- 
ing venting causes mass to be transferred across the liquid -vapor interface. As venting 
continues, it is assumed that only vapor passes through the vent line. The analysis is 
based on a flat interface shape. However, in order to compare the model with the ex- 
perimental data the area of the hemispherical surface is used for the analysis. 

The overall system is divided into two separate control volumes - the vapor and the 
liquid regions. The vapor region consists of the region above the liquid surface, and its 
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volume is a function of the liquid height. The vented vapor M Q is taken from this re- 
gion. The mass transfer terms associated with the vapor region are shown in figure 3. 
As the liquid evaporates during a pressure reduction, some mass M g is transferred to 
the vapor region causing the volume of the liquid region (see fig. 4) to vary slightly with 
time. It is assumed that no boiling occurs within the bulk liquid. 


Derivation of Venting Equations 

In applying the continuity and energy conservation equations, the volume V, which 
appears in the equations, is fixed since actual volume variations due to mass transfer 
will be small. The continuity equation is 



For the vapor region, equation (3) becomes 


dt 




( 2 ) 


(3) 


(4) 


For the liquid region, 


dM^ 

dt 



(5) 


If the internal energy is assumed to be only a function of time, the energy equation in 
integral form is 


- / 
dt Jv 


pu dV + 



puv • n dA = 




n dA 


( 6 ) 
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Using the relation 


h = u + — (7) 

P 

in equation (6) yields the following form of the energy equation with q = 0 (valid only for 
short test times): 

~ f ph dV - — f P dV = - f phv - n dA - f q • n dA (8) 

dtyv dt Jv Ja j a 


For the vapor region, equation (8) becomes 


S'W- V 


dP 


2 _ 


dt 


= M h 


S V, s 


M h 


O V 


(9) 


For the liquid region, 


d 2 

-V, — 2=-M 

dt 1 1 1 dt S C ,*> 


dt 


( 10 ) 


For the bulk liquid region, it is assumed that the initial and final temperatures dur- 
ing venting are identical and equal to Tj. The temperature at the liquid-vapor inter- 
face is Tg sa j., which is the saturation temperature corresponding to the ullage pres- 
sure Pg- The addition of equations (9) and (10), which is not independent of the mass 
transferred across the liquid -vapor interface, yields the product of the interface mass 
transfer rate and the latent heat of vaporization M g hj^. The temperature of the vapor 
is also assumed to be T 2 gat , which is the saturation temperature corresponding to the 
ullage pressure. Since the mass transferred across the liquid-vapor interface enters at 
the saturation temperature, an adiabatic decompression computation was performed to 
show that the vapor initially in the control volume remains nearly saturated during a 
pressure reduction. This analysis was performed for the three fluids employed, and a 
sample calculation is contained in appendix C. The enthalpy of the vapor h y is to be 
evaluated at the saturated conditions. Since the vapor region is of primary interest, 
the governing equations of the vapor region are rewritten as follows: 
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Continuity vapor : 


dM 
v 

dt 



-M 0 


(ID 


Energy conservation vapor : 


£<M^ V ) - = hX - M 0 ) 


( 12 ) 


Expressions for the mass transfer terms which appear in equations (11) and (12) are 
now examined in detail. 

The mass flow rate through the vent M q is determined by using the classical 
choked flow analysis. Since the vented gas goes directly into a vacuum, the choked flow 
'assumption is valid and the mass flow rate is a function of upstream thermodynamic 
properties and an experimentally determined discharge coefficient (see table I). As 
discussed previously, the ullage pressure is a time -dependent quantity. Therefore, it 
fallows that the vent mass flow rate will vary with time. 

i 

Analytical expression for M g . - A conduction analysis by Thomas and Morse 

(ref. 6) was applied to the interfacial mass transfer process occurring during depres- 
surization for an infinitely planer (liquid -vapor) interface. The following expression 
was obtained from this analysis for the mass transfer rate: 





(IS) 


where ^(x) satisfies the transcendental equation 
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TABLE I. - SUMMARY OF PARAMETERS 


Refrig- 

Test 

Initial 

Initial 

Nozzle 

Discharge 

Reduced 

Initial 

Initial 

Final 

Final 

Final 

Dimen- 

Dimen- 

Dimen- 

erant 


filling , 

vapor 

diameter , 

coefficient. 

flov? rate, 

ullage 

ullage 

experi - 

analytical 

adiabatic 

' sionless 

sionless 

sionless 



percent 

volume , 

ni 

C D 

QA. 

pressure, 

temper- 

mental 

ullage 

ullage 

experi- 

analytical 

adiabatic 



liquid 

3 

rn 



1 

p r 

ature, 

ullage 

pressure, 

pressure, 

mental 

pressure 

pressure 







ullage/sec 

N /m 2 

T i- 

pressure, 

N/m 2 

N/m 2 

pressure 

drop, 

drop, 









K 

N/m 2 



drop, 

iP anal /P l 

AP adia/ P l 













AP /P, 
exp’ 1 
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1 

32 

1.93X10' 4 

0. 406X10 “ 3 

0.64 

0.035 

6.96X10 4 

294.3 

8.62X10 4 

8.16X10 4 

7.97X10 4 

0.04 

0.09 

0.11 


2 

29 

2.01 

.889 

.69 

.17 

8.79 

294.7 

7.03 

5.63 

4.94 

.20 

.36 

.44 


3 

33 

1.90 

1.07 

.86 

.33' 

9.10 

293.7 

6.07 

4.07 

3.01 

.33 

.55 

.67 


4 

32 

1.93 

1.32 

.875 

.51 

9.72 

296.5 

5.38 

2.94 

1.63 

.45 

.70 

.83 


5 

32 

1.93 

1,93 

.77 

1.12 

10.1 

295.4 

4.14 

1.31 



.15 

.59 

.88 

.99 

C318 

6 

33 

1. 90x10 " 4 

0. 406x10 “ 3 

0.64 

0.030 

27.9X10 4 

295.9 

26. 9X1 0 4 

25. 5X10 4 

25.25X10 4 

0.04 

0.09 

0.10 


7 

34 

1.87 

.889 

.69 

.16 

30.3 

298.7 

22.1 

19.1 

17.8 

.27 

.37 

.41 


8 

36 

1.81 

.889 

.69 

. 165 

29.0 

297.3 

21.0 

18.0 

16.7 

.28 

.38 

.42 


9 

34 

1.87 

1.07 

.86 

.29 

30.0 

297.3 

17.2 

13.1 

11.4 

.43 

.56 

.62 


10 

35 

1.84 

1.32 

.875 

.455 

29.0 

296.3 

13.0 

8.10 

6.35 

.55 

.72 

.78 

600 

11 

32 

1.93X10" 4 

0. 330X10 “ 3 

0.77 

0.041 

23. 3x10 4 ' 

297.0 

21.7X10 4 

20.8X10 4 

20 . 3X1 0 4 

0.07 

0.11 ' 

0.13 


12 

35 

1.84 

.711 

.81 

.21 

22.8 

294.7 

16.5 

12.8 

11.3 

.28 

.44 

.50 


13 

34 

1.87 

. 889 

.69 

.27 

21.0 

293.7 

14.8 

10.0 

8.50 

.30 

.52 

.60 


14 

35' 

1.84 

.889 

.69 

.28 

24.0 

296.7 

16.4 

11.2 

9.45 

.32 

.53 ' 

.61 


15 

34 

1.87 

1.07 

.86 

.49 

22.8 

296.8 

10.7 

6.30 

4.50 

.53 

.72 

.80 


16 

34 

1.87 

1.93 

.77 

1.31 

23.6 

297.0 

5.65 

1.52 

.63 

.76 

,94 

.97 




where T 1 is the liquid temperature, T 2 t the liquid -vapor interface temperature, 
and Tg y the vapor temperature. Thomas and Morse show that if the expression for 
« z (x) is solved explicitly for ^(x) by expanding both the error functions and exponen- 
tials for small ^(x) one obtains 




K z( T 2,sat- T l) , V T 2,sat- T 2.y) 

V<*1 V«v_ 

V., - 2K v (T 2 ? sat-V> PL + 2K 3 (T 2 >S at-^ 


(15) 


l Pl n fg 




OL 




Now let gat = T 2 ? v ^ as P er our sta ted assumption and appendix C verification): 


2K i (T 2 sat -T 1 ) 


(16) 


Let 


F 1 " V0! lPl h ig 


(17) 


and 


F 2 = 2K Z< T 


2, sat " T l ) 


(18) 


TABLE H. - TYPICAL PROPERTIES FOR TEST FLUIDS AT 297 K 


Refrigerant 

Thermal conductivity 
of liquid, 

K i- 

W/(m)(K) 

Liquid density, 

»v 

kg/m 3 

Thermal diffusivity 
of liquid, 

a l' 

m^/sec 

Heat of vaporization, 

*v 

J/kg 

Ell 

0.0866 

1480 

6. 74X10 ^ 8 

1.84XI0 5 

RC 318 

.0433 

1500 

2.56 

1.07 

R600 

.1056 

575 

7. 88 

3.70 
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Typical values of and Fg and hence the ratio Fg/Fj are now computed. This 
computation was based on the largest temperature change of the tests that were per- 
formed for each fluid. The resulting values of F 2 /Fj were 0.071 for Rll, 0.150 for 
RC318, and 0.157 for R600 (table II). Hence, the largest error in neglecting Fg with 
respect to F^ is 16 percent. Therefore, 

K /( T 2 qat “ T 1 ) 

^(x) a 1 -AS! L (] 

01 1 Pi hfg 

Equation (19) can be substituted into equation (13) to obtain the following approximate 
solution for the interface mass transfer rate : 


M 


s 


A s K !< T l- T 2.sat> 




( 20 ) 


By definition 


K l ~ a l C p, l p l 


( 21 ) 


Therefore, 


M 


s 


A s p Z C p,Z (T l " T 2, sat^ 




( 22 ) 


When combined with the appropriate initial and boundary conditions, equations (11), (12), 
and (22) provide a complete thermodynamic description of the vapor space. These equa- 
tions, along with equation (1) which gives the pressure reduction during an adiabatic de- 
compression, were solved using a computer program (appendix D) which employed a cal- 
culation procedure using 60 time intervals. A comparison of the solutions to these 
equations with experimental data is presented in the RESULTS AND DISCUSSION section. 
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RESULTS AND DISCUSSION 


Pressure Decay Characteristics 

During space missions involving the use of cryogenic liquids, the maximum allow- 
able tank pressure and the thermal environment dictate such items as venting times 
and insulation requirements. Consequently, the pressure response of the tank ullage 
during a venting sequence in space is an important parameter to be considered. A typ- 
ical ullage pressure response for RC318 during one of the experimental tests is shown 
in figure 5. As previously mentioned in the APPARATUS AND PROCEDURE section, 
the experiment tank contents were allowed sufficient time to reach a saturated equilib- 
rium condition corresponding to the ambient temperature in the test facility. The initial 
pressure for the test shown in figure 5 is 29xl0 4 newtons per square meter . A time 
duration of 1.9 seconds was then allowed during free-fall so that the liquid-vapor inter- 
face could attain a hemispherical equilibrium configuration. Note that the tank pressure 
remains constant during this time period as expected. At 1.9 seconds after the initiation 
of the test, the vent was opened and vapor passed through the vent line for a period of 
3 seconds. The final ullage pressure reached 13X10 4 newtons per square meter as seen 
in figure 5. 

In table I, the test parameters for the 16 no-boiling zero-gravity tests conducted 
during this program are shown. Note that test 10 is the typical test whose pressure 
characteristics are shown in figure 5. 


Comparison of Experiment With Theory 

Asa direct result of one of the assumptions made in the analyses, namely, that the 
vapor remains saturated during a depressurization or venting sequence, the mass con- 
servation equation (eq. (11)) and the energy conservation equation (eq. (12)) remain un- 
coupled. Hence, the solution to either equation along with the corresponding expres- 
sions for the vent mass flow rate and interface mass transfer, using the proper bound- 
ary and initial conditions, results in identical solutions. Computer programs were run 
with each separately and yielded similar results. The test conditions and the results 
from the computer program are contained in table I. The final ullage pressure, based 
on the interface mass transfer venting model, is under the column labeled "Final ana- 
lytical ullage pressure. " The tabulated results for the final ullage pressure based on 
the adiabatic decompression venting model are under the column labeled "Final adia- 
batic ullage pressure. " In all cases the adiabatic calculation yielded a lower final 
ullage pressure than for the venting model calculation which included interfacial mass 
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transfer. The difference is due to the fact that the addition of mass from the interface 
into the vapor space would not allow the tank pressure to decrease as rapidly. In all 
tests, for both the adiabatic decompression and the venting model calculation including 
mass transfer, the computed final ullage pressures were lower than the corresponding 
experimental values. 

A graphical comparison of the three previously mentioned pressure drops nondimen- 
sionalized in terms of the initial pressure (AP/P^) is presented as a function of the re- 
duced flow rate (Q/V y ^) in figures 6 to 8. The reduced flow rate is the average volumet- 
ric flow rate divided by the initial vapor volume. (These two parameters are also tab- 
ulated in table I). Figure 6 presents the results for refrigerant 11, figure 7 presents 
the results for refrigerant C318, and figure 8 presents the results for refrigerant 600. 
The experimental and analytical data points contained on these three curves are for the 
0. 06 meter diameter cylindrical tank having a nominal 33 percent by volume initial 
filling of liquid. 

Figures 6 to 8 indicate that the analysis which includes interfacial mass transfer 
shows approximately 25 to 35 percent improvement with respect to the final ullage pres- 
sure prediction in comparison with the adiabatic computation. However, there still 
exists a significant gap between experiment and theory . 

The inclusion of the term F 2 in equation (16) would lead to larger interfacial mass 
transfer rates since the final temperature (T 2 gat ) is always less than the initial tem- 
perature (Tj). Since the discrepancy between experiment and analysis is the result of an 
insufficient amount of interfacial mass transfer, inclusion of the term F 2 would result 
in an improvement. However, this was earlier judged to be a minor effect. Of major 
importance is the fact that an infinitely planer conduction analysis attributable to Thomas 
and Morse (ref. 6) was employed. This, in effect, neglects heat conduction from the 
side walls. Since the liquid -vapor interface was hemispherical in shape there exists a 
very large temperature gradient and, hence, significant heat-transfer rate through the 
thin liquid surface on the walls in the neighborhood of the leading edge. It would appear 
that an improvement in the analysis that included both curved interfaces and finite walls 
would yield higher mass transfer rates and, thus, lead to a much better agreement be- 
tween analysis and experiment. 

Another factor which could contribute to increasing the interfacial mass transfer is 
the liquid motion caused by the initial change from a normal -gravity to a zero-gravity 
configuration and as a response to the venting disturbance. However, it would probably 
be impossible to incorporate this effect into a venting model analysis. 
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Bulk Boiling 


Any complete study of venting under any gravitational conditions should include the 
possibility of boiling occurring within the liquid bulk. The analysis which appears in 
this report does not apply to situations in which bulk boiling occurs. Hence, all the 
data points which appear in table I are for cases wherein no bulk boiling was observed. 
The problem of predicting the inception of bulk boiling is presently without solution. 
Certainly the inception point would depend on the tank surface conditions which determine 
the size and distribution of nucleation sites and also the properties of the particular 
fluid. While the intent of this report was not geared directly toward a discussion of the 
prediction of the aforementioned phenomenon, it did arise in one test case. 

The photographic sequence for test 10 (refrigerant C318) is shown in figure 9. Note 
that the average vent rate for this particular test was Q/V = 0.5 ullage volume per 

second. In this test no bulk boiling occurred. However, as seen in figure 10, for an 
average vent rate of 1.0 ullage volume per second, for RC318, extensive bulk boiling 
occurs. The liquid -vapor interface is pushed toward the vent as a result of the growth 
of two rather large vapor bubbles. These vapor bubbles do not break the liquid -vapor 
interface and cause large surface disturbances as the vent sequence continues. The 
venting model analysis cannot be applied to test cases where bulk boiling occurs since 
the liquid -vapor interface area is not only varying with time but would also be impossible 
to predict. 


CONCLUDING REMARKS 

An experimental investigation of venting cylindrical containers partially filled with 
initially saturated liquids under zero-gravity conditions was conducted in the NASA 
Lewis Research Center 5-second zero-gravity facility. The test fluids, refrigerants 11. 
C318, and 600, all possessed a near -zero contact angle on the container surface resulting 
in a hemispherical liquid-vapor interface shape in zero gravity. The experimental pres- 
sure responses were compared with pressure responses predicted by an adiabatic de- 
compression venting model and a lumped system venting model which includes inter- 
facial mass transfer . 

For the tests which exhibited no bulk boiling the adiabatic decompression venting 
model predicted too large a pressure reduction by nearly a factor of two. The addition 
of the effects of interfacial mass transfer to a lumped system venting model, based on a 
conduction analysis for an infinitely planer (flat) surface, resulted in approximately a 
30 -percent improvement in the determination of the ullage pressure response when com- 
pared to the adiabatic decompression venting model. 
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It is the authors’ belief that the container walls act as a heat source and cause addi> 
tional liquid evaporation thus reducing the experimental pressure decay. In order to 
improve the agreement between experimental and analytical venting system results, the 
infinitely planer surface analytical approach presented in this report would have to be 
modified to account not only for the container walls but also the curvature of the liquid- 
vapor interface. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 14, 1973, 

909-72. 
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APPENDIX A 


APPARATUS AND PROCEDURE 
Test Facility 

The experiment data for this study were obtained in the 5- to 10 -second zero- 
gravity facility at the Lewis Research Center. A schematic diagram of this facility is 
shown in figure 11. The facility consists of a concrete-lined 8. 5 -meter -diameter shaft 
that extends 155 meters below ground level. A steel vacuum chamber, 6.1 meters in 
diameter and 143 meters high, it contained within the concrete shaft. The pressure in 
this vacuum chamber is reduced to 13.3 newtons per square meter by utilizing the 
Center’s wind tunnel exhaust system and an exhauster system located in the facility. 

The ground -level service building has, as its major elements, a shop area, a con- 
trol room, and a clean room. Assembling, servicing, and balancing the experiment 
vehicle are accomplished in the shop area. Tests are conducted from the control room 
(see fig. 12) which contains the exhauster control system, the experiment vehicle pre- 
drop checkout and control system, and die data retrieval system. Those components of 
the experiment which are in contact with the test fluid are prepared in the facility’s class 
10,000 clean room. The major elements of the clean room are an ultrasonic cleaning 
system (fig. 13(a)) and a class 100 laminar-flow station (fig. 13(b)) for preparing those 
experiments requiring more than normal cleanliness. 

Mode of operation . - The zero-gravity facility has two modes of operation. One is 
to allow die experiment vehicle to free-fall from the top of the vacuum chamber, which 
results in nominally 5 seconds of free-fall time. The second mode is to project the ex- 
periment vehicle upwards from the bottom of the vacuum chamber by a high pressure 
pneumatic accelerator located on the vertical axis of the chamber. The total up-and- 
down trajectory of the experiment vehicle results in nominally 10 seconds of free-fall 
time. The 5-second mode of operation was used for this experimental study. 

In either mode of operation, die experiment vehicle falls freely; that is, no guide 
wires, electrical lines, and so forth are connected to the vehicle. Therefore, the only 
force (aside from gravity) acting on the freely falling experiment vehicle is due to re- 
sidual air drag. This results in an equivalent gravitational acceleration acting on the 
experiment which is estimated to be of the order of 10 " 5 g maximum. 

Recovery system . - After the experiment vehicle has traversed the total length of 
the vacuum chamber, it is decelerated in a 3. 6-meter-diameter, 6. 1-meter-deep con- 
tainer which is located on the vertical axis of the chamber and filled with small pellets 
of expanded polystyrene. The deceleration rate (averaging 32 g’s) is controlled by the 
flow of pellets through the area between the experiment vehicle and the wall of the 
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deceleration container. This deceleration container is mounted on a cart which can be 
retracted prior to utilizing the 10 -second mode of operation. In this mode of operation, 
the cart is deployed after the experiment vehicle is projected upward by the pneumatic 
accelerator. The deceleration container mounted on the cart is shown in figure 14. 


Experiment Vehicle 

The' experiment vehicle consisted, of two basic sections (see fig. 15) . The experi- 
ment section is contained in the cylindrical midsection and the telemetry section is con- 
tained in the top fairing. 

Experiment . - The experiment section consisted of the test container tray (see 
fig. 16) plus electrical power and control system equipment mounted in the cylindrical 
section of the experiment vehicle. The test container tray includes the test container, 
camera, and lighting and timing systems. The vent system which included a solenoid 
valve and various sized nozzles and orifices was mounted above the test container. The 
solenoid valve opened during the test drop and vapor was vented either to the low pres- 
sure test chamber or to a collection tank also at a low pressure. The ensuing venting 
procedure was recorded by a high-speed motion picture camera. Elapsed time was ob- 
tained from a digital clock. 

Telemetry system . - The on-board telemetry system is an FM/FM system with 
18 continuous channels. During a test drop, telemetry is used to continuously record 
the output from two low -gravity accelerometers and the tank pressure. The initial tem- 
perature in the vicinity of the liquid -vapor" interface is also obtained from the recording 
system, but the temperature data obtained during the test was not used due to the ques- 
tionable response time of the transducer. The accuracy of all temperatures quoted in 
this report is conservatively estimated to be ±1 . 5°, while the accuracy of all pressures 
is estimated to be ±1 newton per square centimeter. 


Test Procedure 

Cleaning, filling, and hermetic sealing the test containers were conducted in the 
zero gravity facility's clean room (fig. 13). Contamination of the liquid and cylinder, 
which could alter the surface tension and contact angle, was carefully avoided. The test 
cylinders were cleaned ultrasonically in a detergent -water solution, rinsed with a 
distilled -water -methanol solution, and dried in a warm air dryer. The test cylinders 
were rinsed with the test liquid, filled to the desired liquid depth, and sealed to prevent 
contamination. They were then mounted on the test container tray. During the test, a 
predetermined time increment was allowed (1.9 sec maximum) so that the liquid-vapor 
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interface could approach its low -gravity equilibrium shape. After the formation time, 
the solenoid valve was opened and the vapor vented to vacuum for 3 seconds. 

Electrical timers on the experiment vehicle are set to control the initiation and dur- 
ation of all functions programmed during the drop. The experiment vehicle is balanced 
about its vertical axis to ensure an accurate drop trajectory. 

The vehicle is then positioned at the top of the vacuum chamber as shown in fig- 
ure 17. It is suspended by the support shaft on a hinged-plate release mechanism. Dur- 
ing vacuum chamber pumpdown and prior to release, monitoring of experiment vehicle 
systems is accomplished through an umbilical cable attached to the top of the support 
shaft. Electrical power is supplied from ground equipment. The system is then 
switched to internal power a few minutes before release. The umbilical cable is re- 
motely pulled from the support shaft 0. 5 second prior to release. The vehicle is re- 
leased by pneumatically shearing a bolt that holds the hinged plate in the closed position. 
No measurable disturbances are imparted to the experiment vehicle by this release pro- 
cedure. 

The total free-fall test time obtained in this mode of operation is 5. 16 seconds. 
During the test drop, the vehicle’s trajectory and deceleration are monitored in closed- 
circuit television. Following the test drop, the vacuum chamber is vented to the atmos- 
phere and the experiment vehicle is returned to ground level . 
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APPENDIX B 


ADIABATIC DECOMPRESSION VENTING MODEL 

A container of gas at some high pressure is assumed. At time zero, a valve is 
opened allowing the gas to escape. The temperature and pressure dependence with 
time will be examined. An adiabatic process is assumed (q = 0) so that the rate of 
change of internal energy is equal to the rate of energy removal: 


3£.«i(h) 

dt dt 


(Bl) 


where 


U = Mu = MC y T 


(B2) 


and 


h = u + RT 


(B3) 


Expanding equation (B2) gives 


dU 

dt 


= /hl\ + /Vu\ fdM \ 

\at / \dt y \3M/ \dt y 


(B4) 


and 


— = MC v f— ) + u( dM ^ 

dt v \dt / V dt 



(B5) 


When equation (Bl) is used, equation (B5) becomes 


dt v \dt / Vdt; 


(B6) 


The rate of change of mass leaving the container can be expressed in terms of the 
volumetric flow rate Q as follows: 
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dM _ PQ 
dt RT 


when Q is assumed to be a constant in time. Substituting equation (B7) into equa- 
tion (B6) gives 

_ ( u + RT) = MC v f— ^ - — u 
RT V Vdt / RT 


which becomes after some simplification 


Integrating gives 


= __Q_ 

T v v C V 


S. JL 

V v C V 


The solution of the previous equation is 


T 2 = T l e 


-QRt/V v C y 


Now the time dependence of the ullage pressure is derived by employing conservation of 
mass. The rate of decrease of mass in the container equals the withdrawal rate which 
the volumetric removal rate is assumed constant. The rate of change of mass leaving 
the container is expressed in terms of Q by equation (B6) and M, of course, is related 
to the other thermodynamic variables through the perfect gas law 


Therefore, expanding equation (B12) gives 


3P/ \dt 


dM = / 3M \ /dP\ + /9M \ /dT 


3T / Vdt 
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Using previous relations in equation (B12) gives 


PQ _ V v /dP\ _ PV v /dT\ 
RT RT\dt) RT 2\dty 


(B14) 


Rewriting equation (B14) gives 


_§.dt= -4E+ 31 (B16 ) 

v v p T 


Integration yields 


P 1 T 2 

P 2 T 1 


Qt/v_. 


= e 


(B 16 ) 
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APPENDIX C 


EXAMINATION OF ASSUMPTION OF SATURATED VAPOR 

The procedure used to check the validity of the assumption of the vapor temperature 
remaining saturated was to calculate ^ 2^1 from equation (B 11). To accomplish this, 
typical values of Q, V , and t were chosen (table III) . The value of t 2 / t i thus 


TABLE III. - TYPICAL FLUID PROPERTIES AND 


EXPERIMENTAL VALUES FROM TEST 2 


Fluid properties 


Liquid 

Specific heat at constant volume, Cy, J/(kg)(K). 
Gas constant, R, J/(kg)(K) 

Time, t, sec 

3 

Volume of vapor, V y , m 
Average volumetric flow rate, Q, m 3 /sec 
Temperature, T, K 
Initial pressure, P^, N/m 2 

refrigerant 11 
528 
58.6 
3 

2. 01x10 " 4 
0. 342x10 " 4 
294.7 
8.79X10 4 


found was substituted into equation (B16), and equation (B16) was solved for Pg/Pj* 
Finally, the saturation pressure corresponding to Tg is compared with the value of Pg 
computed by means of equation (B16). Using these typical values in equation (Bll) yields 
Tg/Tj = 0.945. Therefore, equals 278. 5 K, which corresponds to a saturation 
pressure of 5.07X10^ newtons per square meter. The ratio P^/P, calculated from 

4*1 

equation (B16) is found to be 0.563. Hence, Pg = 4. 94xl(j newtons per square meter. 
Therefore, the assumption of saturated vapor throughout the venting sequence appears 
reasonable. 
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APPENDIX D 
COMPUTER PROGRAM 

The computer program used for numerically solving equations (11), (12), and (22) 
consists of a main program and four subroutines. The subroutines DD, CURVE, DERIV, 
and RKGS make available the Runge-Kutta method for the numerical solution of an ordi- 
nary first-order differential equation with nonconstant coefficients, the independent 
variable in all cases being time. A calculation procedure was employed which divided 
the 3-second test time into 60 intervals in order to improve the accuracy of the computa- 
tions. Also included in the computer program was the solution of equation (1) which yields 
the pressure decrease based on the adiabatic decompression venting model. This appen- 
dix includes a dictionary of the FORTRAN symbols used, the program listing, and a 
complete input/output list for one test. A typical flow chart is shown in figure 18. All 
of the data were analyzed using the computer program with only slight variations due to 
the different sources and formats of the thermodynamic properties for the three fluids. 


Dictionary of FORTRAN Symbols 


FORTRAN 

Engineering 


symbol 

symbol 


AT 

A 

area of vent nozzle at throat 

AS 

A s 

liquid-vapor surface area 

CD 

C D 

discharge coefficient 

CPL 

C P,1 

specific heat of liquid 

GO 

G o 

acceleration due to gravity 

H 

h 

specific enthalpy 

HFG 

h, 

fg 

heat of vaporization 

K 

K 

thermal conductivity 

MV 

K 

mass flow rate of vented vapor 

MMT 

M s 

mass flow rate across liquid-vapor surface 

P2, Y(l) 

V 

ullage pressure based on venting model 

Q 

Q 

volumetric flow rate 
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FORTRAN 

Engineering 


symbol 

symbol 


T,X 

t 

time 

Tl 

T V,1 

initial ullage temperature 

T2 

T sat, 2 

saturation temperature corresponding to P2 

VV 

V 

volume 

Y3 

^2 

ullage pressure based on adiabatic decompression 

ALPHAL 

a l 

liquid thermal diffusivity 

RHOL 

Pi 

liquid density 

RHOV, Y(2) 

Pv 

vapor density 


PROGRAM LISTING 


EXTERNAL DO 

COMMON P2I »MV,Q, Y3 , AT , A S ,CPL » ALPHAL » VV , RHOL , T1 , CD, GO, KODE 
DIMENSION Y(2) *DY(2) , TITLE (121 
REAL MV 

2 RE AD ( 5 , 102 ) ( TITLE <11 , 1 = 1, 12 I 

READ (5,1051 P2,RH0V,AT,AS,CPL, ALPHAL , VV,RHOL, T1 , CD, GO 
C KODE =1 MMT CALCULATED KOOE=0 MMT = O.D 

K0DE=0 
DO 3 11=1,2 

C INITIALIZE T - X 

X=0»00001 

C INITIALIZE P2 - Y(l) 

Y« 1)=P2 

C INITIALIZE RHOV - Y(2) 

Y ( 2 >=RHOV 
P2I=Y( 1) 

WR I TE ( 6 , 100 ) 

WR ITE< 6,103 ) ITITLEU) ,1=1,121 
WR ITE ( 6, 104) 

WRITE( 6, 107 ) P2, RHOV, AT , AS,CPL, ALPHAL , VV, RHOL , T 1 , CD, GO, KODE 
WRITE (6,104) 

WRITE( 6, 106) 

DO 1 1=1,60 
XF=X+0.05 

CALL RKGS ( X, XF ,H , . IE-3 , Y, DY , 2 , DD) 

WR I TE ( 6, 101 ) X,Y(i),MV,Q,Y3 
1 CONTINUE 
KOOE = KODE-*-l 

3 CONTINUE 
GO TO 2 

IOC FORMAT( 120H1ZER0 GRAVITY VENTING MASS BALANCE AND ENERGY 8ALANCE 
1 LABUS - AMIING ) 

101 FORMAT ( 1CX,5G15. 6) 

102 FORMAT ( 12A6) 

103 FORMAT! 1HK,12A6) 
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onon ono o r> o o o o oooo 


1C4 FORMAT ( 1HK ) 

1C 5 FORMAT (2F6.0»E6.2»7F6.0,F12.0) 

106 FORMAT (lM0»14Xf7HT ( SEC ) » 9X » 2HP2 » 14X , 2HMV » 14X , 1HQ, 14X , 2HY3 1 

107 FORMAT < 5H0P2* ,G16.8/7H RHOV* ,G16.8/5H AT* ,G16 e 8/5H AS* fl G16o8/ 
16H CPl= .G16.8/9H ALPHAL* »G16.8/5H VV= »G16.8/7H RHOL* ,G16.8/5H 
2Tl= * G16. 8/5H CD* ,G16.8/5H GO* ,G16.6/7H KODE* ,12/) 

ENO 


SUBROUTINE OD (X t Y t DY) 

HOMOGENEOUS DIFFERENTIAL EQUATION OF THE FIRST DRDER USING 
RUNGE-KUTTA SOLUTION 

COMMON P2I »MV,Q,Y3,AT, AS, CPL» ALPHAL tW, RHOL, Tl, CD, GO, KODE 
DIMENSION COEFF1 I 5 > ,COEFF2 1 5 > , COEFF3 1 5 ) ,C0EFF4I 5 ) , COEFF5 1 5 > , 
1C0EFF6I5) ,COEFF7<5) ,C0EFFB(5) ,T(100) » P2 A ( 10D > , Y ( 21 ,DY ( 2 ) 

R600 COEFFICIENTS 

CURVE 1 RHOV VS. P 

DATA ICOEFFK I ), 1*1 ,5 ) /0. 1 2 715492E-01 , D. 1D550608E-01 , 0. 11994321 E-0 
14, -C *6138991 IE -06,D. 725379 41 E- 38/ 

CURVE 2 T2 VS. P 

DATA (COEFF21I >, I* l ,5 1 /-45. 905 340, 7. 8234674, -0. 22653 183, 0. 38259717 
IE-02 » -0* 25415846E-04/ 

CURVE 3 H2 VS. P 

DATA ( COEFF3 ( I ) , 1*1 , 5 ) /- 633 . 91 002, 2.41 16720 , -0 • 6 5913 194E-0 1 , 0. 1 550 
134 20E-02.-Q. 6597 1124E-05/ 

CURVE 4 HFG VS. P 

DATA ( C0EFF4I I ) , 1*1 , 5 )/ 1 83. 66624 ,-l. 5806674, 0. 4973961 7E-0 1 , -0.9392 
11366E-03,0.67298545E-05/ 

CURVE 5 K VS, P 

DATA ( COEFF5 ( II , I *i , 5 ) /-0. 5607S 050E-03 , 4. 9051085 , -3. 34892 180E-0 1 , 
10. 719100646-03,-0. 59942675E-05/ 

REAL MMT ,MV,K 
SQRTPI*!. 772438 

COMPUTE MMT - USE THERMODYNAMIC CHARTS TO OBTAIN T2, AND HFG 
P2 = Y( 1 1 

CALL CURVE(COEFF2,P2,T2> 

CALL CURVE (COEFF4,P2, HFG) 

MMT =0. 0 

IF(KODE.EQ.O) GO TO 1 

MMT= AS*RH0L*CPL*(T1-T2> *SQRT I A LPHAL/X ) / < SQRTPI *HFG*60« ) 

1 CONTINUE 

CALL CURVE < C0EFF5 ,P2 , K ) 

MV=CD*AT*K 

COMPUTE CP2/DT USE THERMODYNAMIC CHARTS TO OBTAIN RHOV, RHOVP, AND 
H2P 

CALL CURVE1C0EFFI,P2,RH0V) 

CALL CURVE(C0EFF3,P2,H2) 

CALL DERIV(C0EFF1,P2, RHOVP) 

CALL DERIV(COEFF3,P2,H2P) 

DY(2 ) = UMMT-MV)/VV) 

DY ( 1 ) =DY ( 2 ) /RHOVP 
Q=MV/Y(2) 

Y3*P2I*EXPI-Q*X*1.11 /VV) 

RETURN 

END 
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SUBROUTINE CURV E < COEFF, X » Y ) 

DIMENSION C0EFF(5> 

Y-COEFFIl }+X*(C0EFF( 2 I +X*1 COEFF ( 3 >*X# ( COEFF f 4 ) +X*COEFF ( 5 J I ) ) 

RETURN 

END 


SUBROUTINE DERI V! COEFF »X ,DYDX) 

DIMENSION COEFF! 5 ) 

DYDX=C0EFF!2)*X*( 2.D*COEFF (3)+X*!3.3*COEFF(4)+X*4.0*COEFF(5))) 

RETURN 

END 


PARAMETERS 

XO INITAL VALUE FOR X, RETURNED AS XF NORMALLY 

XF FINAL VALUE FOR X 

HH NORMALLY RETURNED AS STEP SIZE USEO NEAR XF 

RETURNED AS 0 IF MORE THAN 25 HALVINGS 
ERR UPPER ERROR BOUND , USED TO CONTROL STEP SIZE 

Y DEPENDENT VARIABLES, SET TO INITAL VALUES, 

RETURNED AT X=XF NORMALLY 
DERY DERIVATIVES DY/DX 
NDIMI NUMBER OF EQUATIONS IN SYSTEM 

FCT EXTERNAL SUBROUTINE WHICH CALCULATES DERIVATIVES 

CALL IS CALL FCT ( X,Y, DERY ) 

SUBROUTINE RKGS < XO, XF , HH, ERROR, Y, OERY»NO 1M l, FCT I 
DIMENSION Y ( 1 1 , DERY { 1 ) 

DIMENSION AUX(7,20J, A!4) ,B( 4) 

DATA A /. 5, .29289322, 1.7071068, .16666667/ 

DATA B /2.,1.,1.,2./ 

NDIM * NDIMI 
ERR * ERROR 
X = XO 
XEND = XF 

H = ( XEND-X ) #.0625 
CHECK ZERO 
I F ( H.EQ.O. I GO TO 40 
CALL FCT! X , Y, DERY I 


PREPARATIONS OF FIRST RUNGE-KUTTA STEP 
DO 3 1=1, NDIM 
AUX!1,U=Y( l I 
AUX(2,I >*DERYU> 

AUX! 3, I ) =0» 

3 AUX (6, I ) =0. 

H = H + H 
I HLF=-I 
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1 STEPO 
I ENDO 


START OF A RUNGE-KUTTA STEP 

4 IF< (X*H-XEND>*H)7,6»5 

GONE PAST, ADJUST H, SET IEND 

5 HH = H 

H = XENO - X 
IEND = 1 
GO TO T 

ON OR WITHIN TOLERANCE, SET IEND 

6 IEND = 1 

HH * H 

KEEP GOING 

7 I TEST * 0 

9 I STEP= I STEP+1 


START OF INNERMOST RUNGE-KUTTA LOOP 
J-l 

10 A J = A ( J ) 

BJ=B(J) 

CJ = AJ 

IFC J.EQ.4) CJ = .5 
DO II 1=1, NDIM 
R1=H*DERY ( I ) 

R2 = AJ*(Rl-BJ*AUX(6,m 

vm*r(i»+R2 

R2=R2+R2+R2 

11 AUX(6,I)=AUX(6,I>+R2-CJ*Rl 
I FI J-4) 12,15,15 

12 J = J+1 

l FC J~3 ) 13 , 14, 13 

13 X=X*.5*H 

14 CALL FCTCX, Y,DERY> 

GOTO ID 

END OF INNERMOST RUNGE-KUTTA LOOP 


TEST OF ACCURACY 

15 IF( ITE$T)16»16*20 

C IN CASE I TEST=D THERE IS NO POSSIBILITY FOR TESTING OF ACCURACY 

16 DO 17 1=1 , NDIM 

17 AUX{4,n=Yin 
I T E S T = 1 

IST6P=ISTEP+ISTEP-2 
IS I HLF = I HLF + 1 
X = X-H 
H = . 5*H 

DO 19 1=1, NDIM 
VC I > = AUX ( 1, I > 

DERYC 1)=AUX( 2,1 ) 

19 AUX(6,H=AUXI3,n 
GOTO 9 
C 

C IN CASE I TEST= 1 TESTING OF ACCURACY IS POSSIBLE 
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20 I M0D=I STEP/2 

I F C f STEP-IMOD-IMODP 21,23 ,21 

21 CALL FCT 5 X» Y» DERY ) 

DO 22 I = 1 1 MDI M 
AUXI5, H=YU J 

22 AUX 4 T p l 5 =DERY ( I ) 

SOTO 9 

C 

C COMPUTATION OF TEST VALUE DELT 

23 DELT = 0<> 

DO 24 1=1 »NDI M 

24 DELT = DELT ♦ ABS ( AUX C 4, II -Y C 1)1 
moELToLEoERR) GO TO 28 

C 

c ERROR IS TOO GREAT 

I FUHLF.GT.25 ) GO TO 38 
DO 27 1=1 »NDI M 

27 A UX ( 4 » l ) = AUX { 5 « I ) 
ISTEP=ISTEP+ISTEP-4 
X = X-H 

I END=3 
GOTO 18 
C 

C RESULT VALUES ARE GOOD 

28 CALL FCT ( X* Y ? DERY J 
00 29 1=1 fNOIM 
AUXU,n»Y(I > 

AUX<2 e D-DERYU ) 

AUX«3, I 5 = AUX ( 6 9 t S 
YU J =AUX ( 5» I > 

29 D£RY( I ) = AUXt 7» I > 

30 DO 31 I = 1 9 NO I M 
YU ) = A JXt 1, n 

31 DERYU)=AUXC2,n 
IFUEN0J32 ? 32,39 

C 

C INCREMENT GETS DOUBLED 

32 I HLF= I HLF-i 

I STEP= ISTEP/2 
H = H*H 

I F ( IHLFI4, 33*33 

33 I MOD= I STEP/ 2 

I Ft ISTEP-IMOD-IMOD)4,34, 4 

34 I F (OELT«GT . . 02*ERR) GO TO 4 
IHLF = IHLF - 1 

I STEP=I STEP/2 
H = H+H 
GOTO 4 
C 
c 

38 HH = 3. 

GO TO 40 

39 XO * X 

40 RETURN 
END 
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Typical Input/Output 


gravity venting mass balance and energy balance 

VEHICLE w ODEL - ELUIO 9600 t BUT ANE > 


p ?a 

34.8G000C2 

R HDV = 

0.3794COOO 

AT* 

Q.667Q0OCQE-O5 

AS* 

0. 60QC9CC0E-O1 

CPI* 

G, 5 600GCQ0 

alphal 

C » 3CEOOOOOE- 

V v = 

0. 65000GC0E-02 

RH0L = 

35.8DSS999 

Tl* 

76.3G0G0C2 

CD” 

0* 69C00CC0 

GO 3 

25G51. 6GCI 

KODE = 

0 


T <SEC> 

P2 

NV 

Q 

Y3 

0. 500100E-01 

34.29195 

0. 680602 E-03 

0.1819176-02 

34.26353 

0.100010 

33.79093 

0.67U86E-D3 

0 #1019286-02 

33.73536 

0.150010 

33.29687 

0.661693E-Q3 

0.1819386-02 

33.21528 

0.200010 

32.80970 

0.652724E-03 

0.1819466-02 

32.70319 

0.250010 

32.32936 

0.643677E-D3 

0.181954E-02 

32.19896 

0.300010 

31.85576 

0.63475ZE-03 

0 • 181960E-02 

31.70248 

0.350010 

31.38884 

0.625947E-03 

0*1819656-02 

31.21364 

0*40001 0 

30.92852 

0*617261 E-03 

0* 181969E-02 

30.73233 

0.450010 

30.47472 

0.608694E-03 

0. 181973E-02 

30,25843 

0.500010 

30.02737 

0.600243E-33 

0.1819776-02 

29.791B3 

0.550010 

29.58638 

0. 591908E-03 

0. 1819796-02 

29.33243 

0.600010 

29.15167 

0. 583638 E-03 

0.181982E-02 

20.00311 

0.650010 

28.72317 

0. 575589E-93 

0.1819846-02 

20.43476 

0.700010 

28.30080 

D.567584E-03 

0.1019066-02 

27.99620 

0.750010 

27,98447 

0.559698E-03 

0.1019076-02 

27*56456 

0.800010 

27.47410 

0.551921E-03 

0.1619006-02 

27.13950 

0.850010 

27.06962 

0.544251E-03 

0.1019096-02 

26.72100 

0.500010 

26.67094 

0.536&87E-03 

Q.181990E-02 

26.30896 

0.550010 

26.27799 

0.529228E-03 
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Pressure ratio, (P^ - Pf)/P^ 



Time from initiation of test, t p sec 

Figure 5. - Tank pressure response during representative data 
run. Test fluid, refrigerant C318; nozzle diameter, 0.132 
centimeter; average vent rate, 0,5 ullage volume per second. 
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Figure 7. - Effect of volumetric flow rate on pres- 
sure drop for refrigerant C318. 
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Pressure ratio, (Pj - P ^/Pj 
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Reduced flow rate, Q/V Vi 

Figure 6. - Effect of volumetric flow rate on pressure drop 
for refrigerant 11. 
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Figure 8. - Effect of volumetric Now rate on pressure drop 
for refrigerant 600. 




(c) Venting occurring. Time from initiation of test, 4.00 (d) Configuration prior to termination of test. No venting; 

seconds. time from initiation of test, 5.16 seconds. 

Figure 9. - Typical test in which no bulk boiling occurs. Test fluid, refrigerant C318; nozzle diameter, 0.132 centimeter; inital ullage pres- 
sure, 2.9UxI0 5 newtons per square meter; average vent rate, 0.5 ullage volume per second. 
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(a) Formation of zero-gravity equilibrium interface. Venting (b) Boiling occurs near liquid-vapor interface. Time from 
begins; time from initiation of test, 1.90 seconds. initiation of test, 2.70 seconds. 



(c) Vapor generation increases; interface rise toward vent. (d) Configuration prior to termination of test. No venti ng 
Time from initiation of test, 4.00 seconds. occurring; time from initiation of test, 5.17 seconds. 

Figure 10. - Occurrence of bulk boiling in zero-qravity. Test fluid, refrigerant C318, nozzle diameter, 0.193 centimeter; 
initial ullage pressure, 2.76x10^ newtons per square meter; average vent rate, 1.0 ullage volume per second. 
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Figure 11. - Schematic diagram of 5- to 10-second zero-gravity facility. 
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Figure 12. - Control room. 



(a) Ultrasonic cleaning system. 
Figure 13. - Facility clean room. 



(b) Laminar flow work station. 
Figure 13. - Concluded 
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Figure 14. - Deceleration system. 
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Figure 15. - Experiment vehicle. 



Figure 16. - Experiment assembly tray, 
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Figure 17. - Vehicle position prior to release. 



Figure 18. - Computer program flow chart. 
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